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Abstract 

A measurement of elastic deeply virtual Compton scattering 7*p using e^p col- 

lision data recorded with the HI detector at HERA is presented. The analysed data sample 
corresponds to an integrated luminosity of 145 pb^^. The cross section is measured as a 
function of the virtuality of the exchanged photon and the centre-of-mass energy W of 
the 7*p system in the kinematic domain 6.5 < < 80 GeV^, 30 < < 140 GeV and 
|t| < 1 GeV^, where t denotes the squared momentum transfer at the proton vertex. The 
cross section is determined differentially in t for different and W values and exponential 
t-slope parameters are derived. The measurements are compared to a NLO QCD calcula- 
tion based on generalised parton distributions. In the context of the dipole approach, the 
geometric scaling property of the DVCS cross section is studied for different values of t. 
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1 Introduction 



Measurements of inclusive deep-inelastic scattering (DIS) of leptons and nucleons allow the 
extraction of Parton Distribution Functions (PDFs) which describe the fraction of the longitudi- 
nal momentum of the nucleon carried by the quarks, anti-quarks and gluons. A shortfall of this 
approach is that the PDFs contain information neither on the correlations between partons nor 
on their transverse distributions. This missing information can be provided by measurements 
of processes in which the nucleon remains intact and the four momentum transfer squared at 
the nucleon vertex, t, is non-zero [1-6]. The simplest such reaction is deeply virtual Comp- 
ton scattering (DVCS), the diffractive scattering of a virtual photon off a proton 7*p — > 7p. 
In high energy electron-proton collisions at HERA, DVCS is accessed through the reaction 
ep — > e^p [7-9]. This reaction also receives a contribution from the purely electromagnetic 
Bethe-Heitler (BH) process, where the photon is emitted from the electron. The BH cross sec- 
tion is precisely calculable in QED and can be subtracted from the total process rate to extract 
the DVCS cross section. 

Perturbative QCD calculations assume that the DVCS reaction involves two partons in the 
proton which carry different longitudinal and transverse momenta. The difference in longitudi- 
nal momentum of the two involved partons, also called skewing, is a consequence of the mass 
difference between the incoming virtual photon and the outgoing real photon. The skewing can 
be described by introducing generalised parton distributions (GPDs) [1-5], which are functions 
of the two unequal momenta and thus encode information on the longitudinal momentum cor- 
relations of partons. Information on the transverse momentum of partons is incorporated in the 
t-dependence of GPDs [2-5]. The t-dependent functions follow particular equations for their 
evolution as a function of the four momentum transfer squared Q'^ of the exchanged virtual 
photon [3-5]. These evolution equations still need to be tested. 

The DVCS cross section can also be interpreted within the dipole model [10-12]. In this 
picture the virtual photon fluctuates into a colour singlet qq pair (or dipole) of a transverse size 
r ~ which subsequently undergoes hard scattering with the gluons in the proton [13]. 
At very small values of the Bjorken scaling variable x the saturation regime of QCD can be 
reached. In this domain, the gluon density in the proton is so large that non-linear effects 
like gluon recombination tame its growth. In the dipole model approach, the transition to the 
saturation regime is characterised by the so-called saturation scale parametrised here as Qs{x) — 
Qq{xq/x)^^/'^ , where Qq, xq and A are parameters [14]. The transition to saturation occurs 
when Q becomes comparable to Qs{x). An important feature of dipole models that incorporate 
saturation is that the total cross section can be expressed as a function of the single variable r: 

at!M')=<Trj{Tl with T = ^. (1) 

This property, called geometric scaling, has already been observed to hold for the total ep DIS 
cross section [15, 16] as well as in DIS on nuclear targets [17] and in diffractive processes [12]. 
It has also recently been addressed in the context of exclusive processes including DVCS [12] 
and extended to cases with non-zero momentum transfer to the proton [18]. 

This paper presents a new measurement of single and double differential DVCS cross sec- 
tions as a function of Q'^ and the 7*p centre-of-mass energy W . The single differential cross 
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section da/dt is also extracted. The data were recorded in the years 2005 and 2006 with the 
HI detector when HERA collided protons of 920 GeV energy with 27.6 GeV electrons. The 
sample corresponds to an integrated luminosity of 145 pb^\ four times larger than the previ- 
ous HI measurement [8] of DVCS in positron-proton collisions. The measurement is carried 
out in the kinematic range 6.5 < < 80 GeV^ 30 < < 140 GeV and \t\ < 1 GeV^. 
The t-dependence of the DVCS cross section, da/dt, is found to be well approximated by an 
exponential form e^**!*!; this parametrisation is used throughout the paper. The and W de- 
pendences of h are studied. A parametrisation of the observed Q"^ dependence of h is used to 
constrain the normalisation of the pQCD predictions based on GPDs. The validity of the skewed 
evolution equations is tested. The geometric scaling property of DVCS is also investigated and 
the cross section is compared with dipole model predictions. The scaling property is studied for 
the first time for different values of t. 



2 Experimental Conditions and Monte Carlo Simulation 

A detailed description of the HI detector can be found in [19]. Here, only the detector compo- 
nents relevant for the present analysis are described. HI uses a right-handed coordinate system 
with the z axis along the beam direction, the +z or "forward" direction being that of the outgo- 
ing proton beam. The polar angle 9 is defined with respect to the z axis and the pseudo-rapidity 
is given by 77 = —In tan ^^/2. The SpaCal [20], a lead scintillating fibre calorimeter, covers 
the backward region (153° < 9 < 176°). Its energy resolution for electromagnetic showers 
is a{E)/E ~ 7.1%/^E/GeV © 1%. The liquid argon (LAr) calorimeter (4° < ^ < 154°) 
is situated inside a solenoidal magnet. The energy resolution for electromagnetic showers is 
a{E)/E ~ 11%/ ^^/E/GcV as obtained from test beam measurements [21]. The main compo- 
nent of the central tracking detector is the central jet chamber CJC (20° < 9 < 160°) which 
consists of two coaxial cylindrical drift chambers with wires parallel to the beam direction. 
The measurement of charged particle transverse momenta is performed in a magnetic field of 
1.16 T, which is uniform over the full tracker volume. The innermost proportional chamber 
CIP (9° < ^ < 171°) is used in this analysis to complement the CJC in the backward region 
for the reconstruction of the interaction vertex. The forward muon detector (FMD) consists of 
a series of drift chambers covering the range 1.9 < i] < 3.7. Primary particles produced at 
larger r] can be detected indirectly in the FMD if they undergo a secondary scattering with the 
beam pipe or other adjacent material. Therefore, the FMD is used in this analysis to provide 
an additional veto against inelastic or proton dissociative events. The luminosity is determined 
from the rate of Bethe-Heitler processes measured using a calorimeter located close to the beam 
pipe at z = — 103 m in the backward direction. 

A dedicated event trigger was set up for this analysis. It is based on topological and neural 
network algorithms and uses correlations between electromagnetic energy deposits of electrons 
or photons in both the LAr and the SpaCal [26]. The combined trigger efficiency is close to 
100%. 

Monte Carlo (MC) simulations are used to estimate the background contributions and the 
corrections that must be applied to the data to account for the finite acceptance and the resolution 
of the detectors. Elastic DVCS events in ep collisions are generated using the Monte Carlo 
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generator MILOU [22], based on the cross section calculation from [23] and using a t-slope 
parameter b = 5.45 GeV^^, as determined in this analysis (see section [STI ). Inelastic DVCS 
events in which the proton dissociates into a baryonic system Y are also simulated with MILOU 
setting the t-slope bpdiss to 1.2 GeV^^, as determined in a dedicated study (see section[3l). The 
Monte Carlo program COMPTON 2.0 [24] is used to simulate elastic and inelastic BH events. 
The background source of diffractive meson events is simulated using the DIFFVM Monte 
Carlo [25]. All generated events are passed through a detailed simulation of the HI detector 
and are subject to the same reconstruction and analysis chain as the data. 



3 Event Selection 

In elastic DVCS events, the scattered electron and the photon are the only particles that should 
give signals in the detector [8]. The scattered proton escapes undetected through the beam pipe. 
The selection of the analysis event sample requires the scattered electron to be detected in the 
SpaCal and the photon in the LAr. The energy of the scattered electron candidate must be 
greater than 15 GeV The photon is required to have a transverse momentum above 2 GeV 
and a polar angle between 25° and 145°. Events are selected if there are either no tracks at 
all or a single central track which is associated with the scattered electron. In order to reject 
inelastic and proton dissociation events, no further energy deposit in the LAr calorimeter larger 
than 1 GeV is allowed and no activity above the noise level should be present in the FMD. 
The influence of QED radiative corrections is reduced by the requirement that the longitudinal 
momentum balance E — Pz he greater than 45 GeV. Here, E denotes the energy and P^ the 
momentum along the beam axis of all measured final state particles. To enhance the DVCS 
signal with respect to the BH contribution and to ensure a large acceptance, the kinematic 
domain is restricted to 6.5 < < go QeV^ and 30 < W < 140 GeV. 

The selected analysis sample contains 2538 events. It is dominated by elastic DVCS events, 
but also contains contributions from the elastic BH process and from the BH and DVCS pro- 
cesses with proton dissociation, e~p — > e^'jY, where the baryonic system Y of mass My is 
undetected. These background contributions are studied in further detail. A control sample of 
BH events is selected. For this sample, it is required that the electron be detected in the LAr and 
the photon in the SpaCal (see [8] for more details). The COMPTON MC describes accurately 
the normalisation and the shapes of the distributions of the kinematic variables for these events. 
The deviations are within 3%, and this value is used subsequently as an estimate for the sys- 
tematic uncertainty on this contribution. A second control sample dominated by inelastic BH 
and DVCS processes is obtained by selecting events with a signal in the FMD. After subtract- 
ing the inelastic BH contribution, as estimated from the COMPTON MC, this sample allows 
the normalisation of the inelastic DVCS process to be determined. Within the model used in 
MILOU [22], the normalisation of the inelastic contribution is directly related to the exponential 
t-slope parameter. The measured event yield corresponds to an exponential t distribution with 
a slope of 1.2 GeV~^ which is subsequently used in the simulation of inelastic DVCS events. 
The corresponding contribution of proton dissociation in the analysis event sample is found to 
be 16 ± 5%. Other backgrounds from diffractive u and </> production with decay modes to final 
states including photons are estimated to be negligible in the kinematic range of the analysis. 
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Contamination from processes with low multiplicity 7r° production was also investigated and 
found to be negligible. 

The reconstruction method for the kinematic variables Q^, x and W relies on the measured 
polar angles of the final state electron and photon (double angle method) [8]. The variable t is 
approximated by the negative square of the transverse momentum of the outgoing proton. The 
latter is computed from the vector sum of the transverse momenta of the final state photon Pt^ 

and of the scattered electron P^^ : t ~ — iPj^ + PtJ^ . The resolution of the t reconstruction lies 
in the range 0.08 to 0.22 GeV2. 

Distributions of selected kinematic variables are presented in figure [T] for the analysis sam- 
ple. The MC expectations of the different processes are also displayed. Each source is nor- 
malised to the data luminosity. A good description of the shape and normalisation of the mea- 
sured distributions is observed. 



4 Cross Section Determination and Systematic Uncertainties 

The DVCS and BH contributions dominate in the analysis phase space. In addition, an interfer- 
ence term contributes to the cross section due to the identical final states of both processes. In 
the leading twist approximation, the main contribution resulting from the interference of the BH 
and DVCS processes is proportional to the cosine of the azimuthal angle of the photorQ [1,27]. 
Since the present measurement is integrated over this angle, the contribution of the interference 
term is estimated to be small (below 1%). The DVCS cross section, 7*p — > 7p, is evaluated in 
each bin i with the bin centre values Qf,Wi,ti, from the total number A^"*^*^ of data events in the 
analysis sample using the expression 

( jyobs _ jyBH _ jyP-diss^ 
C^DVCsiQi, ti) = ^DVCS ' ^DVCsiQi^ ^i)- (2) 

The other numbers in this equation are calculated using the MC simulations described in sec- 
tion [21 N^^ denotes the number of BH events (elastic and inelastic) reconstructed in bin i, 
jyp-diss number of inelastic DVCS background events, Nf^*^^ the number of DVCS events 
computed from the elastic DVCS MC and cr^^ycs the theoretical DVCS cross section used 
for the generation of DVCS MC events. The measured cross section is thus directly corrected 
for detector inefficiencies and acceptances and is expressed at each bin centre value. 

The mean value of the acceptance, defined as the number of MC events reconstructed in a 
bin divided by the number of events generated in the same bin, is 45% over the whole kinematic 
range and reaches 78% for the highest t bin. The systematic errors of the measured DVCS cross 
section are determined by repeating the analysis after applying to the MC appropriate varia- 
tions for each systematic source. The main contribution comes from the acceptance correction 
factors calculated by varying the t-slope parameter set in the elastic DVCS MC by ±8%. The 
uncertainty on the number of elastic DVCS events lost by the application of the FMD veto is 

'The azimuthal angle of the photon is defined in the proton rest frame as the angle between the plane formed 
by the incoming and scattered electron and that formed by the virtual photon and the scattered proton. 
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modelled by a 4% variation of the FMD efficiency. Both error sources together result in an er- 
ror of 10% on the measured elastic DVCS cross section. The uncertainty related to the inelastic 
DVCS background is estimated from the variation of its t-slope parameter by 25% around the 
nominal value of 6 = 1.2 GeV^^. The resulting error on the cross section amounts to 5% on 
average and reaches 15% at high t. The uncertainties related to trigger efficiency, photon iden- 
tification efficiency, radiative corrections and the subtraction of BH background and luminosity 
measurement are each in the range of 2 to 4%. The total systematic uncertainty of the cross 
section amounts to about 15% and is dominated by correlated errors. 

5 Results and Interpretations 

5.1 Cross Sections and ^-dependence 

The complete DVCS sample is used to extract the W dependence of the DVCS cross section 
expressed at = 8 GeV^ as well as the dependence at = 82 GeV The results are 
displayed in figure |2] and are in agreement within errors with the previous measurements [8,9]. 
The steep rise of the cross section with W is an indication of the presence of a hard underlying 
process [28]. The corresponding cross section measurements are shown in table [T] 

Next, the W dependence of the DVCS cross section is determined for three separate ranges 
of and shown in figure Ha). The corresponding cross section measurements are given in 
table |2l A fit of the form is performed to the cross section in each range. Figure ISj^b) 
presents the 5 values obtained as a function of Q^. It is observed that 5 is independent of 
within the errors. Using the complete analysis sample, the value of 5 expressed at = 8 GeV^ 
is found to be 0.74 ± 0.11 ± 0.16, where the first error is statistical and the second systematic. 

The differential cross section as a function of t is displayed in figures Ufa) and (b) for three 
values of and W, respectively. Fits of the form da/dt ~ e"'''*' are performed taking into 
account the statistical and correlated systematic errors; they describe the data well. The derived 
t-slope parameters b^Q"^) and b{W) are displayed in figures Hfc) and (d), respectively. The 
cross section values and the results for b in each and W bin are given in table |3l This 
analysis extends the study of the evolution of b with Q"^ to larger values than in the previous HI 
measurement [8]. This dependence can be parametrised [30] as 

6(g2) =A{1-B log(gV(2 GeV^))) . (3) 

Fitting this function to the measured b values of the present data and to the value obtained at 
Q2 = 4 GeV2 in the previous HI publication [8] yields A = 6.98 ± 0.54 GeV^ and B = 
0.12 ± 0.03. The systematic errors and their point to point correlations were taken into account 
in the fit, resulting in a correlation coefficient between A and B of pab = 0.92. As shown in 
figure Htc) the fit function provides a good description of the measured b values over the whole 
range. The values of 6 as a function of W are measured for the first time and shown in 
figure Sfd). No significant variation of b with W is observed. 

Using the complete analysis sample, the value of b expressed at = 8 GeV^ is found 
to be 5.45 ± 0.19 ± 0.34 GeV^^, where the first error is statistical and the second system- 
atic. Following [2, 6], this t-slope value can be converted to an average impact parameter of 
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<r\> = 0.65 ± 0.02 fm. It corresponds to the transverse extension of partons, dominated 
by sea quarks and gluons for an average value x = 1.2 10^'^, in the plane perpendicular to the 
direction of motion of the proton. This value is related to the size of the core of the proton with 
no account of the peripheral soft structure. 



5.2 QCD Interpretation in Terms of GPDs 

The determination of b^Q"^) described above can be used to study the Q'^ evolution of the GPDs. 
The DVCS cross section integrated over the momentum transfer t can be written [23] as 

aovcsiQ , W) = TMiQ^) ' 

where lmA{'j*p — > 7p)t=o(<5^, W) is the imaginary part of the 'jp scattering amplitude 

at t = and is a small correction due to the real part of the amplitude. In the following, p 
is determined from dispersion relations [11] to be p = tan(|c<j(Q^)). The coefficient uj{Q'^) 
describes the power governing the W dependence of DVCS at a given Q^. It is taken from the 
corresponding power of the rapid rise of the proton structure function F2 at low x {F2 ~ x^'^) 
[31], assuming that it is sufficiently close to the one in DVCS. In the GPD formalism, the 
amplitude A{YP lp)t=o is directly proportional to the GPDs. As shown in the previous 
section, the Q"^ dependence of the t-slope b is non-negligible. Therefore, the evolution of 
the GPDs themselves is accessed by removing this variation of b{Q'^). For this purpose, the 
dimensionless observable S is defined as 



^Dvcs bm 

Using the parametrisation ^ for b{Q'^), S is then calculated for each Q"^ bin from the cross sec- 
tion measurements of this analysis (table [T]) and from those of the previous HI publication [8]. 
The uncertainties on the parameters A and B of ^ are directly propagated to determine the 
error on b(Q'^) at any given Q'^ value. The results for S are presented in figure Oa) together 
with the prediction of a GPD model [30], based on the PDFs parametrisation given in [32]. It is 
observed that the pQCD skewed evolution equations [3-5] provide a reasonable description of 
the measured weak rise of S with Q'^. 

The magnitude of the skewing effects present in the DVCS process can be extracted by 
constructing the ratio of the imaginary parts of the DVCS and DIS amplitudes. At leading order 
in as, this ratio R = Tm A {YP ~^ 1p)t=o/1^ ^ il*P — ^ 7*p)t=o is equal to the ratio of the 
GPDs to the PDFs. The virtual photon is assumed to be mainly transversely polarised in the 
case of the DVCS process due to the real photon in the final state and therefore has to be taken 
as transversely polarised in the DIS amplitude too. The expression for i? as a function of the 
measured observables can be written as 

^ ^ ^Vnanvcs bjQ^) ^ V^rpycs bjQ^) 

aril* p^X) 7(1+7) aEuFrix, Q^) ^(1 + P^) ' 
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using the relation 0-^(7* p — X) = ATt'^aEMFrix-, Q'^)/Q'^ with a em = 1/137. i? is evaluated 
taking Ft = F2 — Fj^ from the QCD analysis presented in [33] and using the parametrisation ([3]) 
for h{Q'^). The measured values of the ratio R for each Q'^ bin are shown in figure Ob) and 
compared with the calculation based on the GPD model proposed in [30]. The typical values 
of R are around 2, whereas in a model without skewing R would be equal to unity. Therefore, 
the present measurement confirms the large effect of skewing. In GPD models, two different 
effects contribute to skewing [3-5]: the kinematics of the DVCS process and the Q"^ evolution 
of the GPDs. The data are compared to a model which takes only the former effect into account. 
The result of this incomplete model is represented by a dotted line in figure Ob). The present 
measurements show that such an approximation is not sufficient to reproduce the total skewing 
effects observed in the data. 



5.3 Geometric Scaling 

As discussed in sectioniH the dipole model represents another possible theoretical approach to 
describe the DVCS reaction. It is therefore interesting to test if the present DVCS measurements 
obey the geometric scaling laws predicted by such models. 

In the following study parameters of the dipole model are taken from an analysis of the 
total DIS cross section [14, 18]. The saturation scale Qs{x) = Qq{xq/ x)^^^"^ is evaluated using 
Qq = 1 GeV, A = 0.25 and xq = 2.7 10"^. The DVCS cross section measurements listed in 
tableland those from the previous HI publication [8] which are measured at different and 
X = Q'^/W'^ values can be represented as a function of the single variable r (see equation ([T])). 
The result is shown in figure ^a). All of the cross section measurements appear to be well 
aligned on a single curve as a function of r. Therefore the DVCS data are compatible with the 
geometric scaling law. The dipole model of [12, 14] is also represented in figure Oa) and gives 
a good description of the cross section measurements over the complete range of r. 

The dependence of the DVCS cross section on r is also studied at four different values of t. 
For this purpose, the cross section is measured differentially in t for three values of W and two 
ranges of (6.5 < < 11 GeV^ and 11 < < go GeV^), as listed in tableH Keeping 
the same parameters Qo, A and xq as previously defined, each value of the differential cross 
section is again represented as a function of r. The results are shown in figure Ob), together 
with the predictions of the dipole model [12, 14]. For these predictions, the t-dependence is 
factorised out as e^^'*', where the global t-slope parameter b measured in section [5TT] is used. A 
reasonable description of the DVCS cross section values in the four t bins is observed, with the 
same saturation scale Qs{x) used in all cases. 



6 Conclusion 

The cross section for deeply virtual Compton scattering 7*p 7p has been measured with the 
HI detector at HERA. The analysis uses the e~p data recorded in 2005 and 2006 corresponding 
to a luminosity of 145 pb^\ four times larger than in the previous HI publication [8]. The 
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measurement is performed in the kinematic range 6.5 < < 80 GeV^, 30 < W < 140 GeV 
and \t\ < 1 GeV2. 

The W dependence of the cross section is well described by a function W^. No significant 
variation of the exponent 5 as a function of is observed. For the total sample a value 6 = 
0.74 ± 0.11 ± 0.16 is determined. The steep rise of the cross section with W indicates a hard 
underlying process. The t-dependence of the cross section is well described by the form e^^'*' 
with an average slope ofb = 5.45 ± 0.19 ± 0.34 GeV^^. This value corresponds to a transverse 
extension of sea quarks and gluons in the proton of a/< > = 0.65 ± 0.02 fm. The t-slopes 
are determined for the first time differentially in W with no significant dependence observed. 
The study of the dependence of b is extended to significantly larger values compared 
to previous measurements. The slopes found in the present analysis and in the previous HI 
publication are in agreement with a slow decrease of as a function of Q^. 

The measurement of b{Q'^) obtained in the present analysis is used to constrain the normal- 
isation and Q"^ dependence of theoretical predictions based on GPDs. It is found that a GPD 
model reproduces well both the DVCS amplitude and its weak rise with Q^. The skewing ef- 
fects have been investigated and are found to be large, as expected in GPD models. Another 
approach based on a dipole model including saturation effects predicts that the cross section 
can be approximated by a function of the single variable, r = Q^/Ql{x) where Qs{x) is the 
saturation scale. The present measurement of the DVCS cross section is found to be compatible 
with such a geometric scaling using the same parameters as derived from inclusive DIS. For the 
first time, this scaling property is observed for different values of t. 



Acknowledgements 

We are grateful to the HERA machine group whose outstanding efforts have made this ex- 
periment possible. We thank the engineers and technicians for their work in constructing and 
maintaining the HI detector, our funding agencies for financial support, the DESY technical 
staff for continual assistance and the DESY directorate for the hospitality which they extend to 
the non DESY members of the collaboration. We would like to thank Markus Diehl for helpful 
discussions. 



References 

[1] M. Diehl, T. Gousset, B. Pire and J. P Ralston, Phys. Lett. B 411 (1997) 193 
[ |hep-ph/9706 344||. 

[2] M. Burkardt, Int. J. Mod. Phys. A 18 (2003) 173 [hep-ph/0207047). 

[3] M. Diehl, Eur. Phys. J. C 25 (2002) 223 [Erratum-ibid. C 31 (2003) 277] 
[hep-ph /0205208J . 

[4] K. Kumericki, D. Miiller and K. Passek-Kumericki, 'hep-ph/07 03 179 



11 



[5] V. Guzey and T. Teckentrup, Phys. Rev. D 74 (2006) 054027 [hep-ph/0607099|. 

[6] L. Frankfurt, M. Strikman and C. Weiss, Ann. Rev. Nucl. Part. Sci. 55 (2005) 403 
[hep-ph/0507286|. 



[7] C. Adloffet al. [HI Collaboration], Phys. Lett. B 517 (2001) 47 [ |hep-ex/01 070051 1 . 

[8] A. Aktas et al. [HI Collaboration], Eur. Phys. J. C 44 (2005) 1 [hep-ex/0505061]. 

[9] S. Chekanov et al. [ZEUS Collaboration], Phys. Lett. B 573 (2003) 46 [hep-ex/0305028J. 



[10] A. Donnachie and H. G. Dosch, Phys. Lett. B 502 (2001) 74 [hep-ph/0010227| 



[11] L. Favart and M. V. T. Machado, Eur. Phys. J. C 29 (2003) 365 [hep-ph/0302079||. 

[12] C. Marquet and L. Schoeffel, Phys. Lett. B 639 (2006) 471 [hep-ph/0606079J . 

[13] A. H. Mueller, Nucl. Phys. B 335 (1990) 115; 

N. N. Nikolaev and B. G. Zakharov, Z. Phys. C 49 (1991) 607. 



[14] E. lancu, K. Itakura and S. Munier, Phys. Lett. B 590 (2004) 199 [hep-ph/03 10338]. 

[15] K. Golec-Biemat and M. Wusthoff, Phys. Rev. D 59 (1999) 014017 [hep-ph/9807513l. 

[16] A. M. Stasto, K. Golec-Biemat and J. Kwiecihski, Phys. Rev. Lett. 86 (2001) 596 
[hep-ph/0007192]. 

[17] A. Freund, K. Rummukainen, H. Weigert and A. Schafer, Phys. Rev. Lett. 90 (2003) 
222002 [hep-ph/0210139J . 

[18] C. Marquet, R. Peschanski and G. Soyez, hep-ph/0702 17JJ 

[19] I. Abt et al. [HI Collaboration], Nucl. Instrum. Meth. A 386 (1997) 310; 
I. Abt et al. [HI Collaboration], Nucl. Instrum. Meth. A 386 (1997) 348. 

[20] R. D. Appuhn et al. [HI SPACAL Group], Nucl. Instrum. Meth. A 386 (1997) 397. 

[21] B. Andrieu et al. [HI Calorimeter Group], Nucl. Instrum. Meth. A 350 (1994) 57. 

[22] E. Perez, L. Schoeffel and L. Favart, |hep-ph/0411389| 

[23] L. L. Frankfurt, A. Freund and M. Strikman, Phys. Rev. D 58 (1998) 114001 [Erratum- 
ibid. D 59 (1999) 119901] [hep-ph/9710356|. 

[24] A. Courau, S. Kermiche, T. Carli and P. Kessler, Quasi-Real QED Compton Monte Carlo, 
Proceedings of the Workshop on Physics at HERA, Vol. 2: 902-915, Hamburg 1991. 

[25] B. List, A. Mastroberardino, DIFFVM: A Monte Carlo generator for diffractive processes 
in ep scattering, Proceedings of the Monte Carlo Generators for HERA physics, DESY- 
PROC-1999-02,p.396. 



12 



[26] B. Roland, "Mesure de la Diffusion Compton a Haute Virtualite a HERA 11", 
Ph.D. thesis, Universite libre de Bruxelles (2007) (in preparation, to appear at 



http://www-hl.desy.de/psfiles/theses/[ ). 



[27] A. V. Belitsky, D. Miiller and A. Kirchner, Nucl. Phys. B 629 (2002) 323 
[hep-ph/01 12108]. 

[28] I. P Ivanov, N. N. Nikolaev and A. A. Savin, Phys. Part. Nucl. 37 (2006) 1 
[ |hep-ph/0501034| . 

[29] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 



[30] A. Freund, Phys. Rev. D 68 (2003) 096006 [jhep-ph/0306012|. 



[31] C. Mloiietal [HI Collaboration], Phys. Lett. B 520 (2001) 183 [hep-ex/0108035]. 

[32] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P Nadolsky and W. K. Tung, JHEP 0207 
(2002) 012 [hep-ph/020lT95l. 



[33] C. MXoiietal. [HI Collaboration], Eur. Phys. J. C 21 (2001) 33 [hep-ex/0012053|. 



13 



g2 [Gev^; 




(7DVCS [nb] 


8.75 


3.59 ±0.21 ±0.41 


15.5 


1.38 ±0.10 ±0.21 


25 


0.58 ±0.09 ±0.09 


55 


0.13 ±0.03 ±0.04 



W [GeV] 


(7DVCS [nb] 


45 


2.91 ±0.20 ±0.25 


70 


3.96 ±0.32 ±0.37 


90 


4.78 ±0.41 ±0.57 


110 


5.55 ±0.57 ±0.88 


130 


6.56 ±1.17 ±1.77 



Table 1: The DVCS cross section 7*p 'jp, ctdvcs^ a function of for W = 82 GcV and 
as a function of W for Q"^ = 8 GeV^ , both for |t| < 1 GeV^. The first errors are statistical, the 
second systematic. 





ctdvcs [nb] 


W [GeV] 


g2 = 8 GeV^ 


g2 = 15.5 GeV2 


g2 = 25 GeV2 


45 

70 

90 

110 

130 


2.60 ± 0.24 ± 0.24 
3.15 ± 0.40 ± 0.33 
5.25 ± 0.55 ± 0.55 
5.11 ± 0.71 ± 0.76 
5.88 ± 1.89 ± 1.26 


0.94 ±0.10 ±0.10 
1.54 ±0.17 ±0.14 
0.95 ± 0.20 ± 0.17 
1.69 ± 0.31 ± 0.33 
2.06 ± 0.51 ± 0.56 


0.35 ±0.13 ±0.07 
0.36 ±0.10 ±0.05 
0.83 ± 0.18 ± 0.09 
0.90 ± 0.23 ± 0.18 
0.90 ± 0.36 ± 0.32 



Table 2: The DVCS cross section 7*p 7p, aDvcs^ as a function of W for three Q'^ values. 
The first errors are statistical, the second systematic. 
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danvcs/dt nb/GeV^ 


W = 82 GeV 


\t\ [GeV J 


Q2 =8 GeV2 


Q2 =15.5 GeV2 


Q2 =25 GeV2 


0.10 
30 
0.50 
0.80 


13.1 ±1.10 ±1.85 
4.69 ± 0.45 ± 0.55 
1.37 ±0.21 ±0.23 
0.19 ±0.04 ±0.06 


4.37 ±0.47 ±0.86 
1.02 ±0.16 ±0.18 
0.49 ±0.08 ±0.08 
0.12 ±0.02 ±0.02 


1.41 ±0.40 ±0.43 
0.71 ±0.16 ±0.08 
0.28 ±0.07 ±0.04 
0.04 ±0.01 ±0.02 


b [GeV-2] 


5.84 ± 0.30 ± 0.35 


5.16 ±0.26 ±0.30 


5.09 ± 0.55 ± 0.60 




g2 = 10 GeV2 


\t\ [GeV J 


W = A0 GeV 


= 70 GeV 


W = 100 GeV 


0.10 
0.30 
0.50 
0.80 


4.99 ±0.66 ±0.54 
1.45 ±0.29 ±0.18 
0.49 ±0.14 ±0.08 
0.12 ± 0.03 ± 0.03 


7.78 ±0.69 ±0.87 
2.74 ±0.31 ±0.30 
0.81 ±0.14 ±0.11 
0.19 ±0.03 ±0.03 


10.9 ±1.14 ±2.36 
3.47 ± 0.42 ± 0.53 
1.49 ±0.21 ±0.24 
0.19 ± 0.04 ± 0.06 


b [GeV-2] 


5.10 ±0.10 ±0.25 


5.31 ± 0.25 ± 0.27 


5. 18 ± 0.31 ± 0.45 



Table 3: The DVCS cross section 7*p —>■ ■jp, differential in t, daDvcs/dt, for three values of 
at = 82 GeV, and for three values of W atQ^ = 10 GeV^. Results for the corresponding 
t-slope parameters b are given. The first errors are statistical, the second systematic. 





daovcs/dt nb/GeV^ 


g2 = 8 GeV2 


\t\ [GeV'] 


1^ = 40 GeV 


1^ = 70 GeV 


W = 100 GeV 


0.10 
0.30 
0.50 
0.80 


8.10 ±1.22 ±0.82 
2.30 ± 0.54 ± 0.28 
0.45 ±0.22 ±0.10 
0.16 ± 0.06 ± 0.03 


10.0 ±1.30 ±1.27 
4.35 ± 0.63 ± 0.46 
1.08 ±0.27 ±0.17 
0.13 ± 0.06 ± 0.04 


16.0 ±2.11 ±2.74 
5.45 ± 0.80 ± 0.73 
1.96 ±0.41 ±0.35 
0.21 ± 0.09 ± 0.08 




g2 = 20 GeV^ 


t\ [GeV^] 


W = 4:0 GeV 


W = 70 GeV 


W = 100 GeV 


0.10 
0.30 
0.50 
0.80 


1.06 ±0.28 ±0.28 
0.33 ±0.07 ±0.07 
0.22 ±0.06 ±0.06 
0.04 ± 0.01 ± 0.01 


2.38 ±0.29 ±0.26 
0.67 ±0.12 ±0.07 
0.24 ±0.05 ±0.03 
0.07 ± 0.01 ± 0.02 


2.98 ±0.49 ±0.85 
0.89 ±0.17 ±0.17 
0.44 ± 0.08 ± 0.08 
0.06 ± 0.02 ± 0.02 



Table 4: The DVCS cross section j*p •jp, differential in t, daDvcs/dt, for three values of 
W extracted in two Q' intervals: 6.5 < < 11 GeV^ and 11 < < 80 GeV^, corrected to 
the central values of = 8 GeV^ and 20 GeV^, respectively. The first errors are statistical, the 
second systematic. 
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Figure 1 : Distributions of the energy (a) and polar angle (b) of the scattered electron, the energy 
(c) and polar angle (d) of the photon, the electron-photon invariant mass (e) and the proton four 
momentum transfer squared \t\ (f). The data are compared with Monte Carlo expectations for 
elastic DVCS, elastic and inelastic BH and inelastic DVCS (labelled DISS. p). All Monte Carlo 
simulations are normalised according to the luminosity of the data. The open histogram shows 
the total prediction and the shaded band its estimated uncertainty. 
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Figure 2: The DVCS cross section as a function of atW = 82 GeV (a) and as a function of 
W atQ^ = 8 GeV^ (b). The results from the previous HI and ZEUS publications [8,9] based 
on HERA I data are also displayed. The inner error bars represent the statistical errors, the outer 
error bars the statistical and systematic errors added in quadrature. 
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Figure 3: The DVCS cross section as a function of W at three values of Q'^ (a). The solid lines 
represent the results of fits of the form W^. The fitted values of 5{Q'^) are shown in (b). The 
inner error bars represent the statistical errors, the outer error bars the statistical and systematic 
errors added in quadrature. 
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Figure 4: The DVCS cross section, differential in t, for three values of expressed at 
W = 82 GeV (a) and for three values of W at = 10 GeV^ (b). The solid lines in (a) and 
(b) represent the results of fits of the form e"^'*'. The fitted t-slope parameters b(Q'^) are shown 
in (c) together with the t-slope parameters from the previous HI publication [8]. The dashed 
curve in (c) represents the result of a fit to the b{Q'^) values using a phenomenological function 
as described in the text. In (d) the fitted t-slope parameters b(W) are shown. The dashed line in 
(d) corresponds to the average value b = 5.45 GeV^^, obtained from a fit to the complete data 
sample of the present measurement. The inner error bars represent the statistical errors and the 
outer error bars the statistical and systematic errors added in quadrature. 
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Figure 5: The observables S and R (see text), shown as a function of in (a) and (b), re- 
spectively. The results from the previous HI publication [8] based on HERA I data are also 
displayed. The inner error bars represent the statistical errors, the outer error bars the statistical 
and systematic errors added in quadrature. The dashed curves show the predictions of the GPD 
model [30, 32]. In (b), the dotted curve shows the prediction of a GPD model based on an ap- 
proximation where only the kinematical part of the skewing effects are taken into account (see 
text). 
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Figure 6: DVCS cross section measurements as a function of the scaling variable r 
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In (a) the results are shown for the full t range \t\ < 1 GeV^ and in (b) at four 
values of t. The cross section measurements from the previous HI publication [8] are also 
shown in (a). The inner error bars indicate the statistical errors, the outer error bars the statis- 
tical and systematic errors added in quadrature. The dashed curves represent the predictions of 
the dipole model [12, 14]. 
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